Chirality of triangular antiferromagnetic clusters as a qubit 
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We show that the chirahty of triangular antiferromagnetic clusters can be used as a qubit even 
if it is entirely decoupled from the total spin of the cluster. In particular, we estimate the orbital 
moment associated to the chirality, and we show that it can be large enough to allow a direct 
measurement of the chirality with a field perpendicular to the cluster. Consequences for molecular 
magnets are discussed, and an alternative implementation with Cu atoms on a surface is proposed, 
for which one- and two-qubit gates are worked out in detail. Decoherence effects are also discussed. 
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It is widely recognized that the use of quantum prop- 
erties of matter to store and treat information could 
bring exciting new possibilities (see e.g. [l[ and refer- 
ences therein). It is therefore important to find simple 
quantum systems (such as qubits, two- level systems) on 
which some essential basic operations could be made re- 
liably. Requirements include the possibility to prepare 
and measure the system in a specific basis, and also to 
couple subsystems to perform coherent transformations 
of at least two qubits. At the same time, it should be 
free enough from decoherence to allow significant num- 
bers of local operations to be performed, and should also 
be scalable to allow for large arrays of qubits. These 
requirements are very hard to fulfill, and many systems 
lack at least one of them. Solid state implementations 
have the potentiality to be scalable to large size, and 
have been actively studied, using e.g. Josephson junc- 
tions [H, electronic spin in quantum dots 3,Q, nuclear 
spin chains [HI or single- molecule magnets 6|. 

In that respect, antiferromagnetic triangular clusters 
are particularly interesting. First of all, a number of 
experimental realizations are already available. But more 
importantly, the presence of an extra degree of freedom 
in addition to the total spin, the chirality, opens new 
possibilities with respect to purely magnetic qubits. A 
first realization has been explored recently in [7]. There 
it was shown that the chirality induces a spin-electric 
effect whereby the total spin of the cluster can be rotated 
by an external electric field when spin-orbit interaction 
is taken into account. In that scheme, the qubit is still 
essentially magnetic, and the measurement relies on the 
spin magnetic moment of the qubit. 

In this Letter, we discuss the alternative possibility 
of qubits entirely based on the chirality. Getting rid of 
the spin is of course a major advantage regarding de- 
coherence since such a qubit would be less sensitive to 
magnetic noise (see below), but this creates at the same 
time a potential problem for measurement since there is 
no Zeeman coupling of the qubit to an external magnetic 
field any more. As we show, the solution to this problem 



relies on the presence of an orbital moment associated to 
the chirality, as recently pointed out in another context 
in [H . This allows the qubit to be measured as if it had a 
spin, with the important difference however that a mag- 
netic field cannot induce a rotation of the qubit. Note 
that a chiral doublet ground state has been recently re- 
ported in a Dysprosium based molecular magnet [sl . For 
reasons explained below, we will discuss in detail an al- 
ternative implementation with Cu atoms on a surface. 

We note that in the quantum information community, 
the coding of qubits in specific degrees of freedom has 
been investigated in the context of noiseless subsystems 
[iq| . In particular, it has been shown that one can use 
certain degrees of freedom of three spins as a qubit, which 
will be isolated from certain types of noise [11]. In this 
regard, our proposal can be seen as a particularly simple 
realization of a noiseless subsystem, where the protected 
degree of freedom is identified with the physical chirality 
and the logical qubit is directly manipulated and mea- 
sured without any encoding procedure. 




FIG. 1: Schematic view of the two chiral states. The two chi- 
ralities are illustrated by their classical analogs, namely spin 
configurations forming umbrella-like structures with a net mo- 
ment along the x-axis, and with opposite chiralities. They are 
associated with currents circulating in opposite directions. 

The usual starting point to discuss triangular antifer- 
romagnets is the Heisenberg model: 

Hhcis = J1S1.S2 + J2S2'Ss + JsSs.Si — qiibH.S (1) 

where S is the total spin. When Ji = J2 = J3 = J, 
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the ground state is fourfold degenerate, and a convenient 
basis is provided by the simultaneous eigenstates of the 
scalar chirality Si.{S2 x Ss) and of the projection of 
the total spin in an arbitrary direction a: 

\R,cr) = {\ — aaa) -\- uj\a — aa) -\- uj'^\aa — a))/\^ 
|L,cr) = {\ — aaa) -\- uj'^\a — aa) -\- uj\aa — a))/\^ 

where uj = exp(2z7r/3). a = ±1/2 refers to 5" while L 
and R stand for left and right and refer to the chiral- 
ity. To get a twofold qubit, we apply a positive mag- 
netic field H in the direction a, which leaves us with 
a twofold degenerate ground state +1/2), |L, +1/2)} 
separated from the other states by an energy qiibH. The 
two ground states are only distinguished by the chirality 
and thus constitute a non-spin qubit (see Fig.([T]) for a 
pictorial description of these states in the case where the 
magnetic field is in the x direction). This is different 
from the situation where non SU(2) invariant terms in 
the Hamiltonian such as Dzyaloshinskii-Moriya interac- 
tions lift the degeneracy, in which case the two ground 
states differ by the chirality and by the spin 0- 
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FIG. 2: (Color online) Orbital moment associated to the 
chirality as a function of the ratio U/t. Right scale: in 
units of Bohr magneton, assuming the empirical dependence 
t(eV) = 10.67/a(A)2 for Cu 4s orbitals, where a is the inter- 
site distance. Left scale: Opposite of the derivative of the 
energy with the phase calculated at = in units of t. 

Let us first concentrate on the measurement problem. 
Since we have already used the magnetic field to create 
the qubit, it might seem hopeless to rely on the mag- 
netic field to measure it. This is not the case however. 
Indeed, let us turn to a more microscopic description of 
the system in terms of a three-site Hubbard model: 



) + /i.c. + /7^nit^4 (2) 



with implicit periodic boundary conditions. With three 
electrons, a canonical transformation maps this Hamilto- 
nian onto the Heisenberg model with coupling J = At^ /U 



to second order mt/U . However, as recently emphasized 
in [8], charge fiuctuations are not completely suppressed, 
and an orbital current proportional to S'i.(6'2 x S^) can 
be obtained by inverting the canonical transformation. 
These orbital currents create an orbital moment perpen- 
dicular to the plane of the triangle, which can couple to 
an external magnetic field. Let us look at the magni- 
tude of this orbital moment. In [H, this was calculated 
in perturbation theory and it was found that it is of order 
jU^ ^ hence very small. However, having in mind atomic 
Cu, where the outer electron resides in the extended 4s 
orbital, we have calculated the orbital moment for an ar- 
bitrary value of the ratio tjU . The orbital moment is 
related to the magnetic field by /ior6 = —dE/dB\B^Q- 
The magnetic field is taken perpendicular to the plane 
of the triangle and is treated with the Peierls substitu- 
tion as a phase: t = \t\e^'^ . The phase is such that its 
circulation is equal to the fiux of the field, which leads 
to Lp = (it /2\/^)Bo? / where a is the intersite distance 
and 00 = hc/e is the fiux quantum. Using all symmetries, 
the ground state energy can be obtained as: 



E = p cos 
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with p = ^/4^P/9TT2t^, A = -9/p{U^ + 27t^), B = 
2/73/27 and C = eV^t^. For small and large U/t, 
this leads to the following asymptotic expressions for the 
derivative: 



dE 

dip 



\/3t(l - C/V6i^), U/t<€.l 
18^/3^Vf^^ U/t > 1 



(4) 



The full dependence as a function of U/t is depicted 
in Fig.(|2]). The moment is a relatively slowly decreas- 
ing function of U/t: at /7/t = 5, it is still equal to 
about half its value at = 0. To get an estimate of 
the actual magnitude of the moment, we use Harrisson's 
empirical dependence of the hopping integral on the in- 
tersite distance [l2| relevant for the 4s orbitals of Cu: 
t{eV) = 10.67/a(A)^. This leads to an orbital moment 
jiiorb = 0.65/iB at U = 0, a surprisingly large value in 
view of the results of Ref. (sf. For realistic values of /7/t, 
the orbital moment is thus expected to be a sizable frac- 
tion of a Bohr magneton (see Fig.(|2])). So one can use 
a field parallel to the plane of the triangle to create the 
twofold degenerate ground state, and measure the chiral- 
ity by inducing a component of the field perpendicular to 
the cluster to lift the degeneracy. 

This effect will presumably be hard to detect in cur- 
rently available Cus clusters for two reasons: first of all, 
the orbital moment depends on the magnitude of the 
inter-site hopping, which will be much smaller for 3d elec- 
trons than for 45 electrons, as assumed in the above es- 
timate. Besides, if the qubit is of mixed character due to 
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spin-orbit coupling, an external magnetic field will cou- 
ple to both the orbital moment and the total spin, and 
the spin moment will dominate the splitting. 

We thus turn to a description of an alternative imple- 
mentation of such a qubit (see Fig.(j3j) for a representa- 
tion of a chain of such qubits). The basic idea is to work 
with atomic Cu deposited on a surface rather than with 
molecules with Cu^+ ions. The main advantages are: i) 
The absence of spin-orbit coupling since the electron that 
carries the spin is in a 4s orbital; ii) The extension of the 
4s orbitals, which leads to a larger ratio t/U . The impo- 
sition of a magnetic field perpendicular to the plane of 
the triangles will align the qubit towards a state whose 
orbital moment is polarized in the direction of the field. 
This allows to prepare the qubit in a given chirality. Po- 
sitioning of atoms in well-defined arrays with precision 
below the Angstrom level has been routinely realized us- 
ing STM tips to move and place atoms (see e.g. [ill for a 
review) . Equilateral triangular arrays are especially easy 
to build taking advantage of the crystalline structure of 
the atoms of the substrate. 




FIG. 3: (Color online) An example of chirality-based qubit 
chain: triangular clusters of Cu atoms deposited on a surface. 
J is the interaction strength between spins inside one triangle. 
Its spatial modulation induced by approaching an STM tip of 
a triangle controls the one-qubit gates (see text). J' is the 
interaction strength between closest spins of two neighboring 
triangles induced by approaching an STM tip. It controls the 
two-qubit gates (see text). H is the in-plane magnetic field. 

To use this system as a qubit, one should be able to 
manipulate this degree of freedom. To create one-qubit 
gates and modify the state of one qubit, one can modify 
the couplings in a controlled way by use of STM tips (in 
a way which has some similarity with the ones proposed 
in 0, Q ) • The control can be achieved by calibrating the 
system, performing repeated preliminary modifications of 
the couplings and measuring the result. If the couplings 
are transformed from Ji = J2 = J3 = J to Ji = J + JJi, 
J2 = J + ^J2, J3 = J + (^Js, the states \L) and \R) are 
not any more eigenstates of ([1]). The Hamiltonian ([1]) in 
the basis (|L), reads: 



H 



-D - 3 J/4 Vi + iV2 \ 
Vi - iV2 -D - 3 J/4 ) 



(5) 



where the real parameters Vi and V2 are defined by 
D = -{SJi + SJ2 + SJs)/4. and Vi + iV2 = {ujSJi + SJ2 + 
uj'^SJs)/2. The evolution operator through the pertur- 
bation is therefore U = ex.p—i{DJ + Vit^ — ^2^^)^ = 



ex.p{—iDt) exp— i(Vir^ — 1^2^^)^, where / is the identity 
operator and and are Pauli matrices associated 
with the chirality. By an appropriate choice of the three 
parameters (5Ji, SJ2^ SJs and of the time of evolution t, it 
is therefore possible to engineer any rotation of the qubit 
of axis X OTV. It is known that any one-qubit gate can 
be written [l| as ex.p{ia)Ry{f3)Rx{j)Ry{S)^ where Ry{0) 
(resp. RxiO)) is the rotation of angle 9 around the ?/-axis 
(resp. X-axis), and a,/3,7,(5 are real parameters. Thus 
one can generate any one-qubit gate by a sequence of well 
chosen modifications of the couplings inside the triangle. 

To perform a two-qubit gate, one should couple two 
neighboring triangles. If the triangles are sufficiently far 
apart, the permanent coupling between the spins of dif- 
ferent triangles can be made small enough to be negli- 
gible. However, if there are not too far apart, a STM 
tip approached between the closest spins of two neigh- 
boring triangles can mediate an interaction between two 
spins. There is obviously an optimal distance between 
the triangles to respect in the best possible way these 
two requirements. In the presence of the STM tip, the 
new interaction J' between the closest spins of each tri- 
angle creates an interaction between the two triangles i 
and j which, in the basis of the qubits, takes the form 

H' = {J' /^%.Sj{l + 2rf + 2V^rf )(1 - 4r/), (6) 

where Si and Sj are the total spins of the triangles (lij . 
Since the spins are polarized by the magnetic field par- 
allel to the surface, we can consider Si.Sj = 1/4 to be 
constant. Since an arbitrary one-qubit gate can be made 
on each triangle, this can be transformed by rotations of 
the qubits to (J736)(l - 4rf )(1 - 4r/). The evolution 
operator through this interaction is up to one-qubit gates 
exp(-z(4J79)rf r|^t); if applied for a time t = 97r/(16J0 
it generates the two-qubit gate exp(— Z7r/4)diag(l, i, i, 1) 
in the basis (|00), |01), |10), |11)), equivalent up to one- 
qubit operations to a CNOT gate Since the CNOT 
gate plus one-qubit gates form a universal set [l^, the 
model proposed can generate any arbitrary sequence of 
gates of a quantum algorithm. 

The interaction time necessary to implement the one- 
and two-qubit gates is tg ex 1/(5 J, where 5 J is the typ- 
ical magnitude of the process (\uj5Ji + 5J2 + uj'^5J^\ in 
Eq.([5]), 5 J' in Eq.®). If the tip is mechanically moved, 
time scales reached in experiments of the early 90 's were 
of the order of 10 ns [2V\. This time can be signifi- 
cantly reduced, down to the picosecond range, by us- 
ing a non-moving photosensible tip optically addressed 
[22[. The range 10 ps < < 10 ns corresponds to 
0.004 K < ^J < 4 K. With a coupling J of the order 
of 100 K, typical for kinetic exchange, this corresponds 
to small, hence reasonable, modifications of the exchange 
integrals that can be induced by the infiuence of the tip. 
Note that the three important steps in the manipulation 
of a qubit can be optimized independently since they de- 
pend on different parameters: the orbital moment used 
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to initialize and measure the qubit depends on t and U 
in (j2j) and is thus related to J in ([1]), one-qubit gates 
depend on the modulations JJ^, and two-qubit gates are 
controlled by the coupling J' induced between the qubit s. 
This should help to tune the system to the optimal work- 
ing point. We also note that the system operated this way 
should be scalable when increasing the number of qubit s: 
indeed, the qubits are individually addressed and well- 
separated spatially. 

Finally, let us discuss the problem of decoherence and 
why we think such a qubit is well protected. There are a 
priori two sources of decoherence: local vibrations, and 
magnetic noise. Local vibrations can be a source of de- 
coherence because they change the bond lengths inside 
the Cu3 cluster, hence the exchange integrals. As usual, 
zero-point vibrations are not a source of decoherence. In 
the present case, this is easy to see since, as long as the 
cluster has a C3 axis, the ground state for a given total 
spin is two-fold degenerate, even if spin-phonon coupling 
is included. Indeed the degeneracy of the ground state 
on which the qubit is based comes ultimately from the 
fact that the C3 group has a two-dimensional irreducible 
representation to which the ground state belongs, and 
this will remain true for the system including phonons. 

Regarding thermal fluctuations, the main mechanism 
of decoherence comes from the vibrations of individual 
Cu atoms in the potential well in which they are located. 
The depth of the well varies from one system to the other, 
but typical values of the energy barriers for diffusion are 
a fraction of an eV[l8|. The corresponding frequencies of 
harmonic vibrations are in the terahertz range, i.e. 50 K. 
So, by working at low enough temperature, these ther- 
mal vibrations can be exponentially suppressed. Since 
subkelvin temperatures are now accessible to STM ex- 
periments, it should be possible to reduce this source of 
decoherence very efficiently. The other mechanism is re- 
lated to vibrations of the substrate, or phonons. We be- 
lieve that our qubit is naturally protected against these 
vibrations. Indeed, at low temperature, only long- wave 
length acoustic modes can be excited, and they are essen- 
tially decoupled from chirality since it is only sensitive to 
local differences in bond lengths inside the Cus cluster. 

The other source of decoherence is magnetic noise. 
The specificity of the qubit proposed in the present pa- 
per is that the chirality is not coupled to a uniform 
field, as stated before. So decoherence can only be in- 
duced by fluctuating fields that are inhomogeneous on 
the scale of the Cus cluster. A simple calculation simi- 
lar to the derivation of Eq. (|6]) shows that the transition 
amplitude between the two chirality states is equal to 
-(l/3)(/i? + uj'^h^ + ujh^), where a is the direction of 
the field used to fix the orientation of the total spin par- 
allel to the plane, and hi is the local field at site i of 
the cluster. Now, local spins such as nuclear spins in the 
neighborhood of the cluster will produce inhomogeneous 
fields. The typical decoherence time can be expected to 



be similar to that observed in nanomagnets, where such 
a mechanism dominates, hence to be of the order of the 
microsecond [l9|. However, by choosing a substrate with 
little or no nuclear spins, one could get rid of this effect 
altogether. Interesting possibilities are silicon based sub- 
strates such as silicon itself [2o[ or Si02, where the only 
nuclear spins are carried by ^^Si, of natural abundance 
4.7 %, or carbon based substrates such as graphene, with 
nuclear spins on ^^C only, of natural abundance 1 %. A 
few percent of nuclear spins is probably acceptable since 
spins far from the cluster will produce an almost homo- 
geneous field, but in any case isotopic purification could 
be used to reduce further this channel of decoherence. 

So, working at low enough temperature and with the 
appropriate substrate should allow to reach decoherence 
times of several microseconds, much longer than the oper- 
ation time, which is in the nanosecond-picosecond range. 

In conclusion, we have shown how to build a non-spin 
qubit out of a magnetic cluster. The two states of a 
qubit are the lowest energy states of opposite chirality 
of an equilateral triangular cluster of atoms with an in- 
plane magnetic field. We have shown that electric fields 
generated by STM tips are sufficient to create a universal 
set of quantum gates. The presence of an orbital mag- 
netic moment perpendicular to the plane of the triangle 
enables to prepare the system and then to measure it at 
the end. The structure of the chiral states makes them 
robust against noise, and decoherence times should be 
long enough to perform many gates using current tech- 
nology. We thus think our proposal has several attractive 
features for qubit implementation. 
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